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Abstract 

Ratios involving on-resonance measurements of the three-fold correlation 
(TC) and five-fold correlation (FC) cross-sections for which the dependence on 
some of the unknown spectroscopic data is eliminated are considered. Closed- 
form expressions are derived for the statistical distributions of these ratios. 
Implications for bounds on the variance of matrix elements of time-reversal 
(T) non-invariant nucleon-nucleon (NN) interactions are discussed within a 
Bayesian framework and the competitiveness with bounds from other experi- 
ments is evaluated. The prospects for null FC measurements improving by an 
order of magnitude or more upon the current bound on a parity-conserving 
T-odd interaction are good. 
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On-resonance three-fold correlation (TC) and five-fold correlation (FC) neutron transmis- 
sion tests of time-reversal invariance are attractive because of large compound nucleus 
enhancements in sensitivity |5|-|9| . Technically, the experiments have proved challenging be- 
cause of the need to eliminate fake signals and the fact that they require, in addition to a 
polarized neutron beam, polarized and aligned targets, respectively []TU^P~5|] . However, an 



aligned 165 Ho target for FC measurements has been constructed at TUNL 0] and the first 
on-resonance TC measurement with a polarized 139 La target is currently being attempted 
at KEK. 

The complexity of compound nucleus (CN) resonances necessitates a statistical descrip- 
tion of on-resonance correlation cross-sections. The first step towards a quantitative analysis 
of on-resonance FC (TC) measurements is the extraction of the local root-mean-square (rms) 
of CN matrix elements of a parity (P) conserving (non-conserving) T-odd NN interaction 
y(T) (y( PT )) 

or, in the event of null measurements, a bound on this root-mean-square WB 



To this end, it is necessary to introduce a statistical reaction model which relates the un- 
known pure imaginary rms matrix element (V£J2 = ^ v t or HmP = *^pt) to the distribution 
of values of the on-resonance correlation cross-sections. In addition, one has to take into 
account the statistical errors in measurements. Since one is dealing an intrinsically stochas- 
tic observable, the presence of errors can lead to non-trivial restrictions on the minimum 
number M min of on-resonance measurements required in effect to set a bound on and 

VPT 
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To date, consideration of the statistics of on-resonance TC and FC cross-sections has 



been limited to the case of measurements at p-wave resonances [17]]. In the case of the FC 
test, measurements at s-wave resonances are an alternative which can benefit from similar 
enhancements in sensitivity |9j and have the practical advantage that s-wave resonances of 
known spin J have been located in 166 Ho whereas no p-wave resonances have been found 
yet. Recently fTS] , it has been suggested that it may be advantageous to study at weak 



s-wave resonances the ratio -Rfd of (in the notation of [ID]) the FC cross-section o\i2 to 



the deformation effect cross-section (T022 (which is of interest in its own right [20fl) because 



dependence on unknown spectroscopic parameters is partially eliminated. By the same 
token, it is appropriate to consider the ratio i?TP of the TC cross-section cr m to the parity- 
violating cross-section 0101 at p-wave resonances. (For completeness, we note that the cross- 
sections <JkKA are related to elastic elements Spy ij of the scattering matrix by 



CkKA = 27tA" 



,kKA 2 
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where A is the reduced wavelength, k = (2k + l) 1 / 2 and we adopt standard notations for 
the Clebsch-Gordan, Racah and 9-j coefficients. Via the optical theorem and a partial 
wave expansion of the elastic scattering amplitude, the total cross-section cr to t for neutrons 
incident on a target nucleus can expressed as a linear combination of the cross-sections <j kKA 
given in [fRS.) 

In this paper, we discuss the statistics of null measurements of i?FD and -Rtp- The 
dominant contribution to the cross-sections (T122 and er m at a resonance arises from the 
admixture under the corresponding symmetry- violating interaction with the nearest adjacent 
resonance of the right spin and parity and it has been customary in most previous theoretical 
studies to include only this contribution. In fact, this two resonance approximation has 
proved to be very good in the interpretation of parity- violation data. Inclusion of the 
smaller admixtures from other more distant resonances is possible either along the lines of 



the approximate analytical treatment of [17] or via Monte Carlo simulations (also conducted 
in [|T7|) employing the random Hamiltonian matrix ensembles underlying statistical nuclear 
theory, but, in both approaches, the results are rather clumbersome. On the basis of the 



findings of jI7], it can be anticipated that the effect of the admixtures from more distant 
resonances will be to smooth out and broaden the probability density distributions for on- 
resonances values of (T122 and cr m , enhancing, in particular, the tails of the distributions or 
the probability of outliers. For the purpose of the present feasibility study, we choose to 
ignore such effects, making our estimates conservative. 



3 



We show that, in the limit when mixing with only the nearest adjacent resonance is 
taken into account, simple closed-form expressions for the probability densities of -Rfd and 
i?TP can be derived. Motivated by the result of ||16|| , we then investigate the impact of the 
size M of a sample of measurements of -Rfd or -Rtp on the determination of bounds on rms 
matrix elements. Our statistical analysis is performed within a Bayesian framework PI1J22I 
and constitutes the most comprehensive yet of on-resonance reaction tests of time- reversal. 
We end by considering the competitiveness of on-resonance TC and FC measurements with 
other approaches to constraining the T-odd part of the NN interaction. We find that the 
prospects are good for null FC measurements improving by an order of magnitude or more 
upon the bound on extracted recently from neutron-proton charge symmetry breaking 
(CSB) experiments p3| (currently, the tightest bound). 

At a weak s-wave resonance, (T122 is dominated by the contribution arising from the 
mixing under the interaction of the small d-wave component of this resonance (labelled 
1 below) with the s-wave component of the nearest adjacent s-wave resonance of the same 
spin J (labelled 2 below). Retaining only this contribution to cr 12 2 and approximating cr 22 
at resonance 1 by its resonant contribution, the ratio -Rfd = °"i22 / ^022 is given by |TB 



12 /rvT im 



-Rfd = \ =^ z 
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where, in terms of the ratio w = g n ^/ 2 ) / 9 n (s/2) °^ ^he reduced d-wave neutron par- 
tial width amplitudes g^Xx of resonance 1 and the ratio of Racah coefficients p(J) = 
W{J\I2;I\) /W(J\I2-I\) (I is the target spin), 



1 + x l\p{J)w 



1 - y/lp(J)w 

Ek (r n fe) denotes the energy (neutron partial width) of resonance k and V\ 2 denotes the 
pure imaginary matrix element of coupling resonances 1 and 2. Using the known 
spectroscopic data, we work below with tfd = -Rfd/Cfd, which, to the extent that the two- 
level approximation inherent in Eq. (P is valid, possesses the ergodicity property required 
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of observables in statistical nuclear theory f24}] . Unlike a^2, ?"fd (or -Rfd) does not depend 
on the unobservable d-wave neutron partial width of resonance 1. 

The factor z in Eq. ([]]) must be characterized statistically because it involves the ratio of 
reduced partial width amplitudes w which cannot be fixed experimentally. According to the 
empirically grounded Porter-Thomas hypothesis, the reduced amplitudes 9^/2) anc ^ 5n(5/2) 
of (spin J) s-wave resonances are independent gaussian random variables with vanishing 
means and variances given by the running averages (over the resonances) ((SW3/2)) 2 ) an d 
((fin°(5/2)) 2 )' respectively. Thus, the probability density function (pdf) of w is the Cauchy 



distribution p5j [7rA(l + [w/A] 2 )] 1 , where the scale parameter A 



({9. 



(°) )2\ //( a (0) )2 \ 



Using the inverse relation w = w(z) implied by Eq. 

Vi\z) 



we find that the pdf for z 
dw 



dz 



ttA (1 + [w(z)/\] 2 ) 

is also a Cauchy distribution. Below, for simplicity, we use the limiting form 

nr 1 



(3) 
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V 3tt 2 (1 + \z 2 ) 

obtained when / — > 00 and we set A = 1 (taking advantage of the fact that the averages 
((fin°(3/2)) 2 ) an d ((fi , i°(5/2)) 2 ) snomc i be approximately equal). 
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Eq. ® 



Like the CN matrix elements of a generic two-body interaction ||24j| , Im 
can be assumed to be a gaussian random variable of zero mean for different pairs of neigh- 
bouring s-wave resonances of the same J (and with d-wave components). Since the rms 



value of Im 



V (T) 
v \2 



is fx, T FD = z Im 



(T) 
12 



has pdf 



^fd(^fd) = — /(^fd/^t), 



(4) 



where, in terms of the pdf p(z) in Eq. @, 

m -= 1 



^7T 



p I (x/y)e y2/2 



dy 



Substituting for p/, /(x) evaluates to 



foo(x) 



1 1 
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where E\ denotes the exponential integral of order one (defined in Eq. 5.1.1 in p6f). A 
simliar but more complicated expression involving E\ holds for fix) in the general case. We 
have checked numerically for A = 1 that foo(x) does not differ by more than a percent or so 
from the exact result for f(x). 

The dominant contribution to the parity-violating cross-section o"ioi at a p-wave reso- 
nance comes from the mixing of this resonance (labelled p below) via the parity-violating 
interaction V^ p ' with the nearest s-wave resonance of the same spin (labelled s below) ||^| . 
Likewise, the dominant contribution to the TC cross-section am (at the same p-wave res- 
onance) stems from the mixing of these two resonances under the interaction y( PT ) JF|. 
Keeping only these two contributions, the ratio R^p = cr ul /a W i reduces to 

Im 



2/ + 1 . , 
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T , V (P) 
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CTPsin^p + ^Im V<P 



where, in terms of the neutron partial width amplitudes g n (j) an d the total neutron width 
T np of the p-wave resonance, ip p = arctan (g n (3/2) / g n (i/2)) and Xy 2 = g n (i/2) 



arctan 



\/2W{J\ll] ll)/W{J\ll; If) J and V {PT) (V {p) ) denotes the pure imaginary (real) 
matrix element of y( pT ) (V^) coupling the s- and p-wave resonances. 

Since the product Xi^V^ can be extracted (when the two-level approximation holds) 
from the non-zero value of aioi measured at the p-wave resonance (but not x\/2 and 
individually), our candidate for ergodic observable is r TP = i? TP /C T p. The angle (p p for 
different p-wave resonances has (under the Porter-Thomas hypothesis) a uniform random 
distribution, implying that £ = sin(yj p + 5j) has pdf (n^l — C 2 ) 1 - Assuming Im Vi PT ' 
for different pairs of neighbouring s- and p-wave resonances (of the same J) is a gaussian 
random variable of zero mean and rms vpt, ttp = CI m Vpf 7 ^ has pdf 

1 



P Tp (r TP ) 



-g(r T p/v PT ) 



in which, using Eq. 3.383.3 of |2jJ 



j(x) = -^—e-* 2 /*K (x 2 /4), 

7TV Z7T 



where K denotes a zero order modified Bessel function (section 9.6 in p6| ). 

We now consider the implications of our results for the extraction of bounds on 
and vpt using Bayesian methods. To this end, the pdf's Pfd and Ptp are interpreted as 
the conditional pdf's -Pfd^fdI^t) and -Ptp^tpI^pt), respectively. The likelihood function 
for a data set {r{\ of M values of rpp, [r^p] subject to known experimental errors {a{\ is 

M 

L ({ri}\{(Ji} , v) = TT p{ r i\ a h v )i where v denotes [upr] and we assume that the (condi- 

i=l 

tional) pdf p(ri\<Ji,v) for the datum is related to V(r\v) = PpT>(r\v) [Pip(r\v)] by 

+00 

p( ri \<Ti,v) = -f= [ e- {r '- r)2/ ^V(r\v)dr. (5) 



(We adopt the notation and terminology of [21], in which the term likelihood and the 
symbol L are used for both Fisher's likelihood and the probability of a data sample — see 
the footonote on p. 614 of |2l]| .) The corresponding Bayesian posterior for v (= vt,vpt) is 



of the form 

P(v\{r l })=^- 1 L({r l }\{a l },v)p(v), 

where p(v) is the prior distribution for v (which is discussed below) and the constant factor 
M is chosen so that P(v\{ri}) is normalized. 

For the sake of illustration, we deal with data sets where each r{ = 0, which are conceiv- 
ably the best for establishing bounds on vt and vpt- The convolutions in Eq. (||) can then 
be obtained in closed form. Substituting the limiting form of / in Eq. (H), we find (with 
Eqs. 6.225.1 and 6.225.2 in |7|) that, for v T ^ j2/3a tJ 



+00 

/ e- r2 ^P FD (r\v T )dr = - ° l K(^2v T /a l ) } (6) 



n 2 - hi 2 



where K(x) = arccos(x) [arccosh(x)] for x < 1 [x > 1]; for t>x = y2/3<7j, the integral evalu- 
ates to 2/V, the limiting value of the r.h.s. of Eq. @ in this case. The second convolution 
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where we have used Eqs. 6.621.3 and 8.113.1 in but taken the parametric dependence 
of the complete elliptic integral of the first kind K to be given by Eq. 17.3.1 of For 
simplicity, we assume below that all errors take on the same value cr CX p- 

Making use only of the information that vt and vpt are, by definition, non-negative, we 
take the prior p(v) oc 0(f), where G denotes the Heaviside step-function. This choice of 
improper prior is arguably too conservative, but it guarantees that, even for the small data 
samples we anticipate, it is the likelihood L which dominates the posterior P and not the 



prior p (the "jury principle", p. 23 of [p2| ). It is because we are able through our choice 
of prior to rigorously impose the constraint that vt and fpx are non-negative that we have 
adopted a Bayesian approach. 

In Table H, we summarize the information contained in our posteriors P(v\{ri}) by quot- 
ing the bounds v^ a ' on vt and Vpt derived from the relation 

„(«) 



J P{v\{ ri })dv = a 



for three standard choices of the probability a (corresponding to what for a gaussian ran- 
dom variable would be bounds at the 1, 2 and 3 a level, respectively). Because we work 
within a Bayesian framework, the bound v^ a > on (^pt) has the immediate probabilistic 
interpretation that there is a probability a of (^pt) lying m the interval between and 
v( a \ Complementary information on the mean values of Vt and fpx 

oo 

v = vP(v\{ri})dv 



implied by our choice of posteriors is given in Table 0. In the limit of large v, our posterior 
P(v\{ri}) = 0([lnv/v} M ), and so, for M = 1, it is not normalizable and cannot be used to 
infer bounds v^ a \ while, for M — 2, there is no finite mean value v. 

For M > 1, we find via asymptotic analysis that, to leading order in inverse powers of 
M, v T ~ n/(y/6M) and 

y PT ~ 2/ V ttM ; the corresponding asymptotic results for bounds are 
~ ln[l/ (1 — a)]ir/ (\^6M) and, in terms of the inverse of the error function erf (Eq. 7.1.1 



in |26|), ~ 2 erf 1 (a) / \[M. The more rapid decrease of Vt and v!^ with increasing M is 



already appparent for the small values of M in Tables | and [TT[ It is plausible that null values 
of rpD should constrain more than the same number of null values of rxp would constrain 
fPT (assuming that the experimental errors were identical in both instances): Pfd(^) falls 
off less rapidly with increasing x than P^p(x). 

Our findings on the bounds v <Q - ) are consistent with the general (albeit less quantitative) 
assertions of [|T6 1 . In the absence of any prior information, more than one measurement is 



necessary to set a bound. With our choice of prior one needs at least M min = 4 measurements 
to set a 95.5% probability bound comparable in size to the experimental error. It may, 
perhaps, be unwise to attach too much significance to this particular value of M m [ n given 
its dependence on our choice of prior and the fact that weaker bounds can, of course, 
be set with fewer measurements. However, it would seem self-defeating if, when one has 
struggled to reduce experimental errors down to some level, one is not able to set a high 
probability bound of the same order of magnitude. Also, for four (or more) measurements, 
the likelihood function is sufficiently peaked at small v (< cr cxp ) that, for any prior which is 
essentially constant for v < o" exp , much the same bounds will emerge. Use of a prior with this 
behaviour would seem to be natural for an experiment which seeks to improve significantly 
on the current state of knowledge. 

The strength of the interaction V^ T ' is conventionally expressed in terms of the modulus 
~g p of the ratio of the T-odd to T-even p-exchange coupling constants. Bounds on g p at 
the 95% confidence level are 5.8 x 10^ 2 from MeV-neutron FC experiments [28|, 9.3 x 10 -3 



from atomic electric dipole moment (edm) measurements and 6.7 x 10 3 from neutron- 



proton CSB experiments f23| . For the purpose of order of magnitude estimates p9| , |30 



vt (in meV) ~ l0 3 \/WD~g p , where, in the context of the FC measurements under discussion, 
we identify D as the spacing (in eV) between s-wave resonances of the same spin. Thus, 
to improve upon the CSB result, a FC experiment with a 166 Ho target (D = 4.6 eV) will 
need to establish a 95% probability bound on v t of better than 50 meV. This would appear 
to be feasible. Taking Cfd ~ 10 eV -1 (typical values for weak s-wave resonances in 166 Ho 
range from of order unity to of order 10 (18]]), four null measurements of -Rfd with errors 



of 0.1 will yield a 95% probability bound on vt of about 40 meV; for 12 null measurements 
(with errors of 0.1), the corresponding bound is an order of magnitude smaller. Even if 
the null measurements are somewhat less precise and the values of Cfd smaller, it is still 
possible to match the CSB bound with a large enough sample of null measurements. An 
attractive feature of this experiment is that, in the energy interval above neutron threshold 
for which high neutron fluxes are available (up to about 1 keV), there are as many as a 
couple of hundred or so s-wave resonances in 166 Ho which could display a deformation effect 
and hence be used for on-resonance measurements of -Rfd- 

By contrast, it will be significantly more difficult for a TC experiment to reduce the 
already stringent bound on V^ PT \ Neutron and atomic edm measurements imply |31 



an 



upper limit of ~ 10 -5 on the ratio of the coupling strengths of P-odd, T-odd 7r-exchange 
(assumed to be the dominant contribution to y( pT )) to the coupling strength of P-odd, 
T-even p-meson exchange. It is reasonable to assume that this upper limit translates into 
a comparable upper limit on the ratio vp^/vp of rms matrix elements (vp denotes the rms 
matrix element of V^ p '). With Ctp ~ Vp 1 and errors of 1 x 10 -6 in the measurement of 
i?TP (there exist facilities where a statistical error of this magnitude could conceivably be 
attained), four null measurements of i?TP would imply a 95% probability bound of 5 x 10~ 6 
on vpt/vp. 
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TABLES 

TABLE I. Bounds on and vpt in units of the experimental error a c 



M a = ai a a = a.2 a = a% a = a.\ a = a<i a = Q3 

2 7.12 >100 >5000 7.31 >100 >5000 

3 1.53 10.3 80.4 2.09 11.4 81.3 

4 0.786 3.77 17.4 1.30 4.83 19.8 

5 0.521 2.12 7.68 0.995 3.08 9.56 

6 0.387 1.44 4.52 0.828 2.31 6.10 

7 0.308 1.08 3.09 0.721 1.89 4.48 

8 0.255 0.860 2.31 0.645 1.62 3.57 

9 0.218 0.712 1.83 0.589 1.43 2.99 

10 0.190 0.607 1.50 0.545 1.29 2.59 
a «i = 0.683, a 2 = 0.955, a 3 = 0.997. 

TABLE II. Mean values v of vt and vpt in units of the experimental error a cxp . 

M34 5 6 7 8 9 10 

U T 2.89 0.997 0.578 0.403 0.308 0.249 0.209 0.180 

Upt 3.33 1.43 0.968 0.761 0.641 0.563 0.506 0.463 
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